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Recent developments in the CTEQ-TEA global QCD analysis are presented. The parton 
distribution functions CT10-NNLO are described, constructed by comparing data from 
many experiments to NNLO approximations of QCD. 

The global analysis of QCD makes use of experimental data from many short-distance scattering 
processes to construct, within some approximations, universal parton distribution functions 
(PDFs) for the proton. Then these functions can be used to calculate hadronic cross sections in 
the Standard Model and other theories. Global analysis and the resulting PDFs are necessary 
for the interpretation of experimental results at hadron colliders. 

Recently published PDFs are based on next-to-next-to-leading order (NNLO) approxima- 
tions for perturbative QCD [lj. Complete calculations for this order of approximation are avail- 
able for the running coupling ot s (Q), PDF evolution in Q, matrix elements in deep-inelastic 
scattering [2] and vector boson production [3 j . The CTEQ analysis treats quark-mass effects in 
the S-ACOT-x factorization scheme, which has been recently extended to two-loop, or NNLO, 
accuracy [3]. Though the NNLO matrix elements are still unknown for some important pro- 
cesses, such as the inclusive jet production in pp/pp collisions, it is important to use NNLO 
approximations, where available. 

CTEQ has developed PDFs for general-purpose computations and estimates of PDF-driven 
uncertainties over many years |S]. The most recent PDFs in this class, named CT10 and CT10W, 
were published in 2010 [6 . We now present a new family of CTEQ parton distributions, named 
CT10 NNLO. There are several reasons for publishing them. First, the CT10 NNLO global 
analysis is based on the NNLO approximation of perturbative QCD, whereas the CT10 and 
earlier analyses were based on NLO. Second, benchmarking of NLO jet cross sections [7jand 
DIS cross sections was performed to quantify theoretical uncertainties, and an in-depth study 
of the treatment of correlated experimental errors has been completed. Third, selection of 
experimental data sets has been revisited. The new NNLO PDFs are closely related to both 
CT10 and CT10W NLO PDFs and can be matched to either of two NLO PDF sets when 
comparing the NLO and NNLO cross sections. In all three cases, only data from pre-LHC 
experiments were used in the global fit. The same values of the QCD coupling and heavy- 
quark masses as in CT10 NLO were assumed. Some results concerning CT10 NNLO PDFs 
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x f(x,Q) versus x 



were presented at DIS2012 [8_ and will be described here. A longer paper on CT10 NNLO is 
in preparation. The CT10 NNLO PDFs are now available in the LHAPDF library. 

In the near future, a new release of NLO and NNLO 
PDFs. named CT12, will include available data from 
LHC experiments. Some preliminary results of the 
CT12 analysis were also presented at DIS2012 0. 

Selection of data. At NLO, the main distinction 
between CT10 and CT10W sets concerns the inclu- 
sion of the DO Run-2 W electron and muon asymmetry 
Ag(yi), that constrains the behavior of the ratio 
d{x,Q/u(x,Q) at x > 0.1. The CT10 NLO set does 
not include the DO Run-2 A e data, while the CT10W 
NLO set includes 4 p T e bins of A e . The CT10 NNLO 
analysis includes all data sets that were used in the 
NLO fits, with the exception of the Tevatron Run-1 in- 
clusive jet cross sections [9] that have been superceded 
by more precise Tevatron Run-2 jet cross sections |10| ; 
and the DO Run-2 Ag data sets, of which only most 
inclusive (best understood) bins of pti are included 
in both the electron and muon channel. Since CT10 
NNLO includes only a part of the DO Ai data that 
distinguishes between CT10 NLO and CT10W NLO, 
it can be treated as a counterpart of either the CT10 
NLO or CT10W NLO PDF set. 

Overview of the PDFs. Figure [l] gives an 
overview of the CT10 NNLO PDFs. Four PDFs are 
shown: u V2l \ cncc (x , Q) = (u-u)(x,Q); d va ,i ence {x, Q) = 
(d — d)(x,Q); g(x,Q): and q se a{x, Q) = 2(d + u + 
s)(x,Q). The vertical axis is xf(x,Q). The CT10 
NNLO PDFs are illustrated by plotting all the error 
PDFs; hence the figure shows not only the central fit 
but also the uncertainty ranges. The dashed curves are 
the central-fit CT10 NLO PDFs. 

Both NLO and NNLO fits have about the same 
X 2 /N pt w 1.1 for N pt = 2700 data points. Slide 5 
in Ref. [8] shows a more complete comparison of CT10 
NNLO to CT10W NLO, for Q = 2GeV and for thre 
parton flavors, g, u, and u. The various PDFs are plot- 
ted as a ratio to the central CT10W NLO. The curves are the ratios of the central CT10 NNLO 
to CT10W NLO. The shaded regions are the error bands for the PDFs (both NLO and NNLO). 
The central NNLO PDFs differ from the central NLO PDFs, but the difference is comparable 
in size to the error bands. The error band for NNLO is slightly smaller than for NLO. 

Compared to CT10W NLO, the NNLO PDF set at a small scale Q has a suppressed gluon 
and increased sea quarks at x < 10 -2 , reduced g(x, Q) and d(x, Q) at x > 0.1, and very different 
charm and bottom PDFs (slide 6 in [8 ). The reduction in g(x,Q) Compared to MSTW'08 
NNLO, the central CT10 NNLO gluon PDF is somewhat harder at x < 10~ 3 and x = 0.1 - 0.5, 
and softer at x > 0.5 (slide 9 in [8]). The strangeness PDF is larger at x ~ 10~ 2 in CT10 
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Figure 1: CT10 NNLO (solid color) 
and NLO (dashed) parton distribution 
functions. 
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NNLO than in MSTW'08 NLO, producing a good agreement with the ATLAS measurement of 
the s(x)/u(x) at this x value. 



Boson/collider 


CT10 NLO 


CT10 NNLO 


MSTW'08 NNLO 


W+ LHC14 (nb) 


12.2 ±0.5 


12.7 ±0.5 


12.4 ±0.2 


W+ LHC7 (nb) 


6.0 ±0.2 


6.3 ±0.2 


6.2 ±0.1 


W + Tevatron (nb) 


1.35 ±0.05 


1.38 ±0.05 


1.38 ±0.02 


W~ LHCT4 (nb) 


8.9 ±0.4 


9.4 ±0.4 


9.3 ±0.2 


W~ LHC'7 (nb) 


4.10 ±0.15 


4.29 ±0.16 


4.31 ± 0.07 


Z LHC14 (nb) 


2.07 ±0.08 


2. 17 ±0.08 


2.13 ±0.03 


Z LHC7 (nb) 


0.96 ±0.03 


1.00 ±0.03 


0.99 ±0.02 


Z Tevatron (pb) 


260 ±9 


263 ±8 


261 ±5 


Hg M LHC14 (pb) 


101 ±9 


99 ±8 


102 ±7 


H" SM LHC7 (pb) 


31.2 ± 1.9 


29.7 ± 1.7 


29.8 ± 1.3 


Hg M Tevatron (pb) 


1.77 ±0.12 


1.77 ±0.12 


1.80 ± 0.11 



Table 1: Total cross sections for production of electroweak bosons. 




Figure 2: Comparisons of ATLAS data with ResBos predictions for Z° and W lepton rapidity 
distributions. 

Predictions for the LHC. In a future paper we will provide detailed comparisons of 
theory and data, where the theory is calculated from the CT10 NNLO PDFs. Here we collect 
some representative cross sections for the hadron colliders. Table [T] compares predictions for 
total cross sections for W, Z and Higgs boson production via gluon fusion (with Higgs mass of 
125 GeV) at the Tevatron and the LHC (with y/s =7 and 14 TeV). The comparison is between 
CT10 NLO, CT10 NNLO, and MSTW'08 NNLO. The CT10 NNLO central PDF increases 
the total cross sections by a few percent compared with CT10 NLO accuracy and is close to 
MSTW'08. Theoretical uncertainties from alternative PDF sets for CT10 NNLO are similar to 
those for CT10, and in W/Z production they are about twice as those for MSTW'08. 

Fig.|2]shows the comparison of ATLAS data [H] with ResBos [T5] predictions for Z and W- 
lepton rapidity distributions at the LHC (y/s = 7 TeV) using CT10 NNLO PDFs. Theoretical 
uncertainty bands were calculated using the error PDF sets. The ResBos prediction of Z and 
VF + -lepton rapidity distribution, using the central PDF set, is higher than ATLAS data by 
a few percent. However, for VF _ -lepton rapidity distribution, the ResBos prediction is more 
consistent with ATLAS data. It is expected that these data could further refine the PDFs at 
the NNLO accuracy. 
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Figure 3: Comparison of ATLAS data for inclu- 
sive jet pt distribution with a theoretical pre- 
diction using CT10 NNLO. 



Fig. [3] compares the ATLAS data for inclusive jet transverse momentum distribution with 
theoretical predictions based on the NLO matrix elements and CT10 NNLO PDFs. They agree 
well even without including the systematic shifts, except for the large rapidity region. After 
accounting for the systematic shifts, the reduced \ 2 is 0-78 for the measurement with R=0.4 
and 0.76 for the one with R=0.6. The effect of the LHC data on the PDFs will be explored in 
the CT12 analysis. 

This work was supported by the U.S. DOE Early Career Research Award DE-SC0003870 
and by the U.S. NSF under grant No. PHY-0855561. 
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